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Nuclear magnetic resonance (NMR), electron paramagnetic resonance (EPR), magnetization mea- 
surements and electronic structure calculations in VOM0O4 are presented. It is found that VOM0O4 
is a frustrated two-dimensional antiferromagnet on a square lattice with competing exchange inter- 
actions along the side (Ji) and the diagonal (J2) of the square. From magnetization measurements 
Ji + J2 is estimated around 155 K, in satisfactory agreement with the values derived from electronic 
structure calculations. Around 100 K a structural distortion, possibly driven by the frustration, is 
evidenced. This distortion induces significant modifications in the NMR and EPR spectra which can 
be accounted for by valence fluctuations. The analysis of the spectra suggests that the size of the 
domains where the lattice is distorted progressively grows on cooling as the temperature approaches 
the transition to the magnetic ground state at Tc — 42 K. 

PACS numbers: 76.60.Es, 76.75.-fi, 75.10.Jm 

I. INTRODUCTION 

In the last decade transition metal oxides have attracted a lot of interest in view of the rich phenomenology induced 
by the strong electronic correlations. The properties of these oxides are rather peculiar once the interaction of the 
electrons with the lattice becomes relevant. Jhis is the driving mechanism of several phenomena as, for example, 
superconductivity, colossal magnetoresistivityEl and Spin-Peierls transitioncl. Recently, the importance of the coupling 
between the electron spin and the lattice have emerged for a new class of materials, the frustrated antiferromagnetsB. 
In this case, the magnetoelastic coupling tends to relieve the degeneracy of the ground-state caused by the frustration of 
the magnetic exchange couplingsQ . This is the situation observed, for example, in Li2VOSi04EfQ which is a frustrated 
two-dimensional S — 1/2 antiferromagnet (2DFQIIAF) on a square lattice with competing exchange interactions 
along the side (Ji) and the diagonal (J2) of the square. This compound has Ji ~ J2 ( Ji + J2 — 8.5 K) and represents 
a prototype of the two-dimensional Ji — J2 model, which was extensively studied from a theoretical point of view in 
the last decadeEl. In the absence of any spin-lattice coupling the ground-state is double degenerate, the two states 
corresponding to coUinear phases (hereafter called I and II) which differ in the orientation of the magnetic wave- 
vectorEl. The jjoagnetoelastic coupling leads to a lattice distortion at Tdist — (Ji + >^2)/2, which affects ^^Si and ^Li 
NMR spectraHQ, and Li2VOSi04 is observed to collapse alwiys in one of thf; two possible ground-states. 

Another 2DFQIIAF, nearly isostructural to Li2VOSi04El, is VOM0O4EJ. The structure of these compounds is 
formed by pyling up layers of SiVOs for the former and of M0VO5 for the latter. These layers contain VO5 pyramids 
separated by (Si,Mo)04 tetrahedra (see Fig. |l|). The only difference is that in Li2VOSi04 a plane 0^^!^ ions is present 
between the SiVOs layers. VOM0O4 has been recently investigated by Shitwaki and coworkersO, which, however, 
have considered it as a protoype of a weakly one-dimensional antiferromagnetll3 instead of a 2DFQHAF, as it will be 
shown in the following. The interest for VOM0O4 stems from the fact that although the structure is very similar to 
the one of Li2VOSi04 the exchange couplings are more than an order of magnitude larger. Thus the comparison of 
the properties of the two systems would allow to understand if frustration is indeed the driving mechanism for the 
observed structural distortions. 

In this manuscript nuclear magnetic resonance (NMR) , electron paramagnetic resonance (EPR) and magnetization 
measurements in VOM0O4 powders are presented. The temperature dependence of the susceptibility evidences that 
also for VOM0O4 Ji ~ J2. At temperatures below Ji + J2 , namely at Tdist — 0.64(Ji -I- J2), VOM0O4 shows a 
lattice distortion possibly driven by the magnetic frustration. The distortion seems to induce valence fluctuations, 
not observed in Li2VOSi04 , with a charge transfer from V''"'" to Mo^+. Moreover, it is found that domains of 
distorted lattice, with a size which progressively grows on cooling, are formed below Tdist- The magnitude of the 
superexchange couplings were estimated from electronic structure calculations and the two-dimensional character of 
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VOM0O4 evidenced. Finally, the role of Mo d orbitals in determining the differences with respect to Li2VOSi04 is 
emphasized. 

The paper is organized as follows: in Sect. II the technical aspects and the experimental results will be shown, while 
in Sect. Ill the analysis of the data, including the electronic structure calculation, the analysis of NMR relaxation 
rates and of the lattice distortion will be presented. The final conclusions are summarized in Sect. IV. 

II. EXPERIMENTAL ASPECTS AND EXPERIMENTAL RESULTS 
A. Sample Preparation, EPR and Magnetization Measurements 

VOM0O4 powders were obtained by solid state reaction starting from a stoichiometric mixture of M0O3 (Aldrich, 
99.5+%), V2O5 (Aldrich, 99.99 %) and V2O3 heated in a vacuum-sealed quartz tube at 675 C for 24 hours. V2O3 
itself was prepared by reducing V2O5 (Aldrich, 99.6+ %) under hydrogen at 800 C. The sample purity was analyzed 
by means of X-ray powder diffraction and all diffraction peaks corresponded to the ones of VOM0O4 (JCPDS file : 
18-1454). Single crystals were prepared by chemical transport reaction starting from a stoichiometric mixture of the 
starting materials and TeC14 (10% in weight). The mixture was sealed under vacuum and heated for 24 hours at 575 
C, then slowly cooled at 10 C/hour down to room temperature. 

EPR spectra were recorded with an X-band spectrometer equipped with a standard microwave cavity and a variable 
temperature device. The measurements were performed both on powders and on a single crystal of volume V < 0.03 
mm"^. The EPR powder spectra are characterized by a lineshape which becomes progressively more asymmetric as the 
temperature is lowered below 130 K (see Fig. ^ and the inset to Fig. ^). These spectra can be quite well simulated 
by considering a temperature independent cylindrical g-tensor with components gc — 1.960 and gab = 1.932. These 
values are consistent with the estimates of g in the framework of the prvstal field approximation, by assuming 
a spin-orbit coupling A = 150 cm^^ and some covalency between V and OlJ. It should be noticed that gc and gab 
values are reversed with respect to what one might expect just by looking at the spectrum at T = 70 K in Fig. |^a. 
In fact, one would be tempted to associate the low- field most intense peak with in grains where i? + c, while 
the less intense one with those grains with H \\ c. However, if this assignment is made the data cannot be fitted 
adequately. The increase in the intensity of the low-field peak with respect to the high-field one, below 130 K (see 
the inset to Fig. ||), has rather to be associated with a faster decrease of the linewidth AH for H \\ c than for + c 
(see Fig. These results can be suitably compared to the ones derived from EPR measurements on a single crystal. 
The crystal was mounted on a sample holder that allowed to rotate the field in the ac plane. V^^ EPR spectra on 
the crystal confirmed that the ^-tensor is practically temperature independent down to Tc ~ 42 K with gc and gab 
identical to the ones derived from the EPR powder spectra. Moreover, the temperature dependence of AH is the 
same found for the powders, characterized first by a decrease on cooling, then a minimum around 60 K, and finally an 
increase as the temperature approaches Tc (see Fig. ||). The temperature dependence of the area of the EPR powder 
spectra, which in principle is proportional to the static uniform susceptibility %, is shown in Fig. ^. Above 100 K the 
temperature dependence is very similar to the one derived for x from magnetization measurements (see paragraph 
below), however, below this temperature a rapid decrease of the EPR intensity is evident (see Fig. ^d), down to Tc- 
Then, below Tc a small signal with a different g, possibly arising from impurities, starts to be detected. 

Magnetization (Af) measurements were performed on VOM0O4 powders using a commercial Quantum Design 
MPMS-XL7 SQUID magnetometer. The temperature dependence of the susceptibility, defined as x = M/H with 
H the intensity of the applied magnetic field, is shown in Fig. ||a. One observes a high temperature Curie- Weiss 
behaviour, a broad maximum around 100 K typical of low-dimensional antiferromagnets and a kink at Tc ~ 42 K, 
which indicates the presence of a phase transition. This trend is the same already observed by Shiozaki et al. (Ref. 
p+l ). At temperatures well above the maximum the susceptibility is given by 

x(T)^j^+Xvv, (1) 

where C is Curie constant, 6 the Curie- Weiss temperature and xvv Van-Vleck susceptibility. In order to estimate 
0, which for a 2DFQIIAF on a square lattice is equal to Ji + J2, one has to determine first the value of xvv- Since 
Xvv does not contribute to the EPR signal its value can be directly determined by plotting x measured with the 
SQUID against the EPR area for T > 150 K (see Fig. |^). One finds xvv = 3.5 x 10"'' emu/mole, a value consistent 
with the separation between the t2g levels derived from crystal field calculations and close to the one estimated for 
Li2VOSi04 , where V*^ has practically the same coordination. Then, by fitting the susceptibility data for T > 150 
K, with Eq. (0) one derives 6 = 155 ± 20 K. 
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B. NMR spectra and relcLxation rates 

^^Mo NMR spectra were recorded both on unoricnted as well as on magnetically aligned powders by summing the 
Fourier transform of half of the echo signal recorded at different frequencies. The powders were oriented in epoxy 
resin with the magnetic field direction along the c axis which, as can be seen from a close inspection of the crystal 
symmetries, corresponds to the principal axis of the electric field gradient (EFG) at ^^Mo nuclei. The NMR spectra 
in the oriented powders are characterized by five well-defined peaks (see Eig. |6p) separated by vq ~ 106 kHz. For 
a cylindrical EFG tensor, as the one of ^^Mo in VOM0O4 , one has thata j/q = SeVzzQi^ — 7oo)/20ft., where Vzz 
is the principal component of the EFG tensor, Q = —0.019 barn ^^Mo electric quadrupole moment and (1 — 700) 
Stcrnhcimer antishielding factor. One can compare the experimental value of i^q with the one derived from an estimate 
of the EFG based on lattice sums, within a point charge approximation, and taking (1 — 700) — 24.861, as estimated 
theoreticallytj. One obtains vq ~ 102.5 kHz, in remarkable agreement with the experimental finding. On the other 
hand, the powder spectra are characterized by a sharp central peak, corresponding to the 1/2 —1/2 transition, and 
by an underlying broad powder spectrum. Below Tc ^'"^Mo NMR powder spectrum broadens, as expected in presence 
of a magnetic order yielding a local field at the nuclei which is randomly oriented with respect to the external field. 

It is interesting to analyze the temperature dependence of the shift of the central line in the oriented powders, 
for H II c, and for the unoriented powders, which probe mainly the shift for i? _L c. As reported in Fig. ^ one 
observes that the resonance frequency of the peak in the unoriented powders has a temperature dependence which is 
exactly the opposite of the one observed for the susceptibility (see Fig. Since the quadrupolar corrections to the 
central line shift are negligible one can assert that the opposite behaviour of these two quantities is due to a negative 
hyperfine coupling constant {A) between ^^Mo nucleus and the 4 nearest neighbour V^^ ions. In fact, the shift of the 
NMR line can be written as 

AK=^^+S, (2) 

with fiB the Bohr magneton and 5 the chemical shift. Hence, by plotting AK vs. x (see Fig. ^) one can derive the 
hyperfine coupling constant A which for T < 100 K turns out Apow — —9 kOe. In Fig. ^ one clearly observes that 
around 100 K there is a sizeable change of slope which has to be associated with a marked increase in the hyperfine 
coupling and suggests that around 100 K significant modifications in the local structure around Mo^+ are taking 
place. Above 110 K one has Apow — —2.5 kOe The shift measurements in the oriented powders yield quantitatively 
similar results (see Fig. |^b), pointing out that the hyperfine coupling is quite isotropic. The values for the component 
of the hyperfine coupling tensor for H \\ c turn out, Ac ~ -11.5 kOe for T < 100 K and Ac ~ -2.7 kOc for T > 110 
K. Finally, it must be mentioned that while above 100 K the width of the central line is temperature independent, 
a sizeable broadening is observed below 100 K, the linewidth for -ff || c increasing from about 2.3 kHz at 106 K to 
about 6.5 kHz at 50 K. This fact suggests an increasing inhomogeneity at the microscopic level. 

Nuclear spin-lattice relaxation rate 1/Ti was measured on the central ^^Mo NMR line by means of a saturation 
recovery pulse sequence. The recovery law was found multiexponential, as expected. Now, the point is whether 
the relaxation process is driven by fluctuations of the hyperfine field or of the EFG. As we shall see later on in the 
discussion of the experimental results, the magnetic relaxation mechanism is the dominant one (Sect. IIIB). Then, 
the recovery law for the nuclear magnetization {t) is 
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The values of 1/Ti derived from the fit of the experimental data with Eq. (||) are reported in Fig. |8[ One notices 
a decrease of the relaxation rate on cooling down to about 90 K, then a plateau and a peak at Tc, as expected for a 
second order phase transition. 

The decay of the echo signal after a 7r/2 — t — tt pulse sequence was observed to be practically exponential. The 
decay of the amplitude of ^^Mo echo signal arises in principle from three different contributions, namely 



E{2t) = £;(o) 



D{2t) X 



X e 



(4) 



The first term D{2t), is the decay associated with ^^Mo nuclear dipole-dipole interaction. The second moment M2 of 
the corresponding frequency distribution was determined on the basis of lattice sums taking into account the natural 
abundance of ^^MoEil. It was found that = 90 s ^ . If one takes this value and assumes a gaussian decay, sizeable 
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deviations from the exponential behavior should be detected, mainly at low temperature. The absence of any evidence 
for such gaussian deviation could stem from the low natural abundance of ^^Mo nuclei which leads,, as for diluted 
nuclear spins, to a dipolar contribution to the echo decay which is neither gaussian nor exponentialEll. Then one has 
to approximate the nuclear dipole contribution to the echo decay with a moment expansion. D(2t) was expanded up 
to the fourth moment and its expression was used to fit the data up to 2t ~ l/{2\/M2). It must be remarked that, 
in anyway, above 120 K the contribution of D{2t) to the echo decay is small with respect to the one due to the third 
term in Eq. (^. 

The second term in Eq. (jj) is Redfield contribution to the echo decay. When just the ceptral transition of an 
1 — 5/2 nucleus is irradiated and in the case of an isotropic spin-lattice relaxation rate one hasEl = g/Ti. The 

third term is the dominant one and originates from a low-frequency dynamics characterized by a correlation time 
which modulates the resonance frequency of the nuclei by Aw'. The observation that the echo decay is exponential 
implies that Alj'tc ^ 1, so that / — 2t/tc and the last term of Eq. (||) becomes ea;p(— 2T/T2), witht3 I/T2 = Aw'^Tc. 
The temperature dependence of I/T2 derived by fitting the echo decay with Eq. (^) is shown in Fig. ^ 

Finally, it must be mentioned that a '^^V NMR signal was detected, with a temperature dependence of 1/Ti and 
of I/T2 very similar to the ones reported in Ref. |l6[ However, the values of the relaxation rates are too small to be 
ascribed to sites in VOM0O4 . The comparison of (1/T2)^ with the second moment derived for ^^V dipole-dipole 
interaction shows that this signal must be due to a few percent of ^^V nuclei, possibly belonging to V^^ impurities. 

III. ANALYSIS OF THE DATA AND DISCUSSION 
A. Electronic structure and superexchange couplings 

The superexchange couplings Ji and J2 were estimated both theoretically, starting from electronic structure calcu- 
lations, as well as experimentally from the temperature dependence of the susceptibility. The electronic structure of 
VOM0O4 was calculated by using the density £iMictional theory (DFT) in the local density approximation (LDA). The 
tight binding linear muffin tin orbital methpdila (LMT047 Stuttgart code) was adopted together with the exchange- 
correlation potential of Perdew and ZungerEj while the lattice parameters were taken from Ref. 

In Fig. |l^ the electronic structure of VOM0O4 and the corresponding density of states (DOS), derived with the 
linear tetrahedron method, are shown. The DOS was checked to have already converged with a mesh of about 858 
irreducible k-points. One can notice that only two relatively narrow bands, well separated from all the others, cross 
the Fermi level (ep), which was set to zero energy. The density of states shows a pronounced feature around ep, i.e. 
in correspondence of these two bands. This feature is more evident in the lower part of Fig. |l^ where these two bands 
and the density of states (DOS) in the energy window (-0.8,0.4) eV are reported. In order to minimize the linearization 
error in this energy window, we placed the linearization energies close to the Fermi level. The eigenvectors of the two 
conduction bands have mainly V dxy character, mixed with some O2 Px/y and, in the case of the lower energy band, 
with some Mo dxy character. At the F point the two conduction bands are, respectively, the bonding (lower energy) 
and antibonding (higher energy) V dxy bands. 

The LDA bands can be understood from a few band tight binding model, as shown in Appendix, and the dispersion 
curve of the two conduction bands can be written in terms of the nearest neighbors (NN) (ti) and next nearest 
neighbors (NNN) (12) hoppings within the [001] plane and of the hopping between adjacent planes {t±) 

e(k) = eo + 2t2[cos(fc:i,a) -I- cos(fcya)] 

+4ii[cos(fca;a/2)cos(fcya/2)] -I- 2t_LC0s(fczc), (5) 

Their values can be estimated from a least square fitting of the calculated band structure. The results are shown 
in Tab. |. The NN hopping, ti, and the NNN hopping, t2, have two contributions of opposite sign. The first one 
originates from the hopping between V and NN O orbitals while the second one from hoppings involving V and NN 
O and Mo orbitals. The sign is different because while the energy of Op orbitals lies below the Fermi level the one of 
Mod orbitals is above (see Appendix) . The contribution coming from Mo depends mainly on the energy of Mo da;^ 
effective orbital, whose energy is affected by the hopping with V (1^2 and, therefore, depends on V-Mo distance. Thus, 
the closer is V to Mo the higher is the energy of Mo dxy and the the smaller the contribution of Mo to t2 and <i. In 
addition, the hopping between Op and Mo dxy states tends to enhance the ratio ti/i2 (see Appendix). 

The hopping integrals can now be used to estimate the exchange couplings among V"'+ spins. VOM0O4 is a half 
filled band Hubbard insulator and, in the limit of strong Coulomb repulsion, the exchange couplings can be expressed 
as Ji = 4:tf/{U — Vi). Here ti are the NN and NNN hoppings, U the on-site Coulosib repulsion and Vi the inter-site 
Coulomb repulsion, which is supposed to be much smaller than U. It was shownEj that typical values of U for the 
vanadates are [/ ~ 4 — 5 eV. So, by taking U '-^ 5eV and neglecting Vi the coupling values shown in Tab. | were 
derived. 
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ti 


t2 


t± 


U 


Ji 


J2 


J± 


J1/J2 


lower band 


-110 


-52 


-1 


5 


110 


22 


< 10^2 


4.5 


higher band 


+135 


-42 


-2 


5 


154 


16 


< 10^2 


10 



TABLE I: Hopping integrals (in meV) for VOM0O4. The Coulomb repulsion is in eV and the exchange coupling costants are 
in K. 



First one notices that J± is at least four orders of magnitude smaller than the in-plane coupling constaatSj pointing 
out that VOM0O4 is a 2D system and not a one-dimensional one as claimed by Shiozaki and coworkersOii3. Second, 
it should be observed that the value of Ji + J2 ranges between 132 and 170 K, in good agreement with the value 
Q = Ji + J2 ^ 155 K derived experimentally for the Curie- Weiss temperature. On the other hand, for both bands 
we find Ji > J2, with J1/J2 around 4.5 for the lower energy band and 10 for the higher energy band. This result, 
however, seems in contrast with the experimental findings. In fact, the temperature dependence of the susceptibility 
and in particular the ratio between Q and the temperature of the maximum in the susceptibility are very similar to 
the ones of Li2VOSi04 , pointing out that also in VOM0O4 J2/J1 — 1- The similarity becomes evident once x is 
plotted as a function of T/9 (see Fig. |ll|). It is interesting to observe that an analogous discrepancy between the 
ratiOrJp/Ji derived experimentally and the one estimated from electronic structure calculation was found by Roesner 
et al.Eil for Li2VOSi04 . Also in that case the estimate of Ji -I- J2 was in good agreement with the experimental one, 
while the value of J2 was found about a factor 10 larger than Ji, the opposite of what happens for VOM0O4 . To 
assure that the estimate of J2 / Ji was not influenced by the method adopted for calculating the band structure , the 
coupling constants were calculated also for Li2VOSi04 and a ratio J2/J1 — 10 was found, in good agreement with 
Roesner et al.Eil results. The big difference in the ratio calculated for Li2VOSi04 and VOM0O4 cannot be associated 
with a difference in the V-0 distances, which are quite similar in both compounds, or with the small rotation of the 
basis of the VO5 pyramid. This difference should rather be ascribed to the role of Mo dxy orbitals in VOM0O4 and 
of Li s orbitals in Li2VOSi04. As already mentioned the hopping between Op- Mo dxy tends to enhance the ratio 
Ji/J2- On the other hand, in Li2VOSi04 the hopping between O p and NN Li s orbitals gives a contribution to ti 
only, which has a sign opposite to the one due to the Yd-Op hopping. Hence the hopping through Li s orbitals reduces 
ti and the ratio Ji/J2- 

Therefore, the observed discrepancies between the experimental and calculated values of J2/J1 cannot originate 
from the method adopted to calculate the band structure but must have another origin. They should rather be 
associated with the simplified expression used to derive the superexchange couplings, where just the on-site repulsion 
U was considered. 



B. Mo relaxation rates and EPR linewidth 

As shown in the previous section ^^Mo echo decay, which probes the very low- frequency dynamics, is characterized 
by two regimes: a high temperature one where I/T2 = Aw'^Tc decreases on cooling and a low temperature one where 
l/r2 increases on approaching Tc from above. This means that the correlation time Tc which describes the dynamics 
decreases on cooling from room temperature down to T ~ 100 K. Which could be the origin of these dynamics? 
One possibility is that ^^Mo echo decay above 100 K is driven by the relaxation of unlike spins, namely, taking into 
account the natural abundance and the magnitude of the nuclear magnetic moments present in VOM0O4 , of ^^V 
spins. Then, Tc = Ti of ^^V and Auj' corresponds to the nuclear dipole coupling between ^^Mo and ^^V nuclei which, 
from lattice sums—fcuHis out Aw' ~ 1510 s"^. Now one can directly estimate ^^V 1/Ti from ^^Mo I/T2 experimental 
data (see Fig. ^BB. 

One observes^^V l/Ti increasing exponentially on decreasing temperature, as one would expect for a correlated 
2DQHAFE3. In fact, the nuclear spin-lattice relaxation rate can be written as 

with 1^(5)1^ the form factor, which gives the hyperfine coupling of the nuclei with the spin excitations at wave- vector 
q, and Saa{q,^L) the component of dynamical structure factor at nuclear Larmor frequency. If scaling arguments 
apply one can express Saa{q,^L) in terms of the in-plane correlation length and, provided that VOM0O4 is in the 
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renormalized classical regime and V hyperfine coupling is mainly on-site, one finds 



1 



(T) cx ^(T) ~ 0.49 X exp{2TTps/T) 



1 T 



(7) 



with ps the spin-stifi^ness. If '''^V 1/Ti data are fitted with this simple expression a poor fitting is obtained. The point 
is that most of the data obtained for ^^V 1/Ti lie in a temperature range where T > Ji + J2 ~ 155 K and scaling 
arguments can no longer be applied. A more accurate quantitative aiialysis can be performed for T 3> Ji + J2, where 
V*"*" spins are uncorrelated. In this temperature limit one can writecil 



725(5+1)^2 



Ti 



(8) 



with ^^V hyperfine coupling constants and oje — \J J^^b /K)\/2zS{S + l)/3 the Heisenberg exchange fre- 
quency, where z = 4 is the number of coupled through Ji or through J2 , to a reference V*^ ion. If one takes 
l/Ti ~ 6-Mis~^ for T ^ Ji + J2 (see Fig. |l^), one derives A±^ ~ 80 kGauss. This is a typical value for V"'+ hyperfine 
couplingcS, supporting the assumption that vanadium nuclear spin-lattice relaxation is driving ^^Mo echo decay. 

The increase in ^^Mo l/r2 on approaching Tc must have a different origin since ^^V 1/Ti is expected to continue 
increasing on cooling and finally diverge at the transition temperature. The change of behaviour around 100 K could 
be ascribed to the onset of a very low- frequency dynamics, which is possibly associated with the motions of domain 



walls separating coUinear I and II domains, as recently observed in Li2VOSi04ti3. 

It is interesting to compare the temperature dependence of ^^V and ^^Mo nuclear spin-lattice relaxation rates (see 
Figs. 8 and 12). One observes that while the former increases on cooling the latter decreases. One could then be 
tempted to associate ^^Mo relaxation to another mechanism, for example a quadrupolar one, where the relaxation 
is due to phononaij. However, the nuclear spin-lattice relaxation rate due to phonons turns out to be an order of 
magnitude smaller tloaji the one derived experimentally if the recovery laws appropriate for a quadrupolar relaxation 
mechanism are usedEZl. On the other hand, if one estimates the value expected for 1/Ti in the assumption of a 
relaxation mechanism driven by V^^ dynamics for T ^ Ji + J2 (Eq. @)), the calculated value turns out to be 
slightly larger than the experimental one. So, it is possible that the in-plane spin correlation causes a decrease of 
^^Mo 1/Ti. This is what is expected if ^^Mo form factor filters out, at least partially, the spin fiuctuations at the 
critical wave- vector. In fact ^^Mo form factor, |A(g)p = [2A{cos{qxa/2) + cos{qya/2))]'^ , is peaked at {qx = 0,qy = 0), 
zero at (7r/a,7r/a) and reaches a reduced value at (7r/a,0) (or (0, 7r/a)), which corresponds to the critical wave-vector 
of the envisaged coUinear ground-state. This situation is very similar to the one found in CFTD, a non-frustrated 
2DQHAFEa. In this system have a form factor similar to the one of ^^Mo in VOM0O4 and 1/Ti was also observed 
first to decrease on cooling for T < J and then to increase. An accurate calculation of the temperature dependence 
of ^^Mo 1/Ti in VOM0O4 goes beyond the aim of this work, since it would require the precise knowledge of the 
temperature dependence of the hyperfine coupling constants between 220 K and Tc- 

It is also instructive to compare ^^Mo l/Ti with the EPR linewidth AH (see Fig. The similarity in the 

temperature dependence of both quantities is striking. Although it is not straightforward to establish a relationship 
between these two quantities, the formcp,probing the spectral density of the 2 spins correlation function while the 
latter of the 4 spins correlation functioned, the physical origin of their behavior is the same. In fact, also the initial 
decrease of AH on cooling has to be associated with the loss of weight of the q ^ diffusive modes and to an 
increase in the spectral weight at (ts/a, 0) (or (0,7r/a)), which finally gives rise, in view of the slowing down of the 
critical fiuctuations, to a peak at TfCJ. As pointed out by Richards and Salamon (Ref. ^9|) the transfer of spectral 
weight from q ~ to the critical wave vector causes also a modification in the angular dependence of AH , with first 
a decrease of AHab/ AH^ and then an increase, exactly as it was found for VOM0O4 (see the inset to Fig. 3). For 
T <C Ji + J2, the EPR linewidth should, scale with the in- plane correlation length and if the same scaling laws derived 
for two-dimensional antiferromagnet£j are used, one should find AH cx Then, by fitting the few experimental 
data in Fig. ^ for T < 55 K and assuming the temperature dependence of ^ given by Kcl (|^) one derives a value for 
2t: Ps around 60 K, well below Ji + J2, as expected for a frustrated 2D antiferromagneiE3. 



C. Frustration driven structural distortion 



The analysis of NMR spectra points out that a local structural distortion around ^^Mo nuclei takes place at 
Tdist — 100 K, yielding a sizeable change of the magnetic hyperfine coupling (see Sect. IIB). The occurrence of a 
structural distortion in VOM0O4 is a natural consequence of the frustration which, in the absence of spin-lattice 
interaction, for J2/J1 — 1 would lead to a double, degenerate ground state down to a temperature where an Ising 
transition to one of the two ground states occurs!^. The effect of the lattice is to relieve the degeneracy among the 
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two ground states, so that the frustrated system always collapses in one of the two states. This is somewhat analogous 
to the Jahn- Teller distortion which relieves the degeneracy among the electronic Iffpis split by the crystal field and 
for this reason some authors have called this distortion the " Spin- Teller" disJtartionuH. Recently, evidence for such a 
distortion in a three dimensional pyrochlore antiferromagnet was presentedH On the other hand, some connection 
with the Spin-Pcicrls distortion is also present. In fact, quite recently Becca and Milaa have shown that for a 2D 
Ji — J2 system the magnetoelastic coupling would induce a distortion which, depending on the value of J2/J1 and 
on their dependence on the lattice parameters, could break either just the rotational invariance or, as for a standard 
Spin-Peierls transition, also the translational invariance. Recent NMR measurements suggest that the breakdown 
of the rotational invariance takes place in Li2VOSi04ta In VOM0O4 , however, the situation is somewhat more 
complicated than the one described by Becca and MilatJ since Ji and J2 show a subtle dependence on Mo position 
and, therefore, cannot be expressed in a simple form in terms of — distance. 

Now, based on simple order of magnitude estimates one can show that the the origin of the lattice distortion in 
Li2VOSi04 and VOM0O4 is the same. In fact, if one takes the ratio between the temperature at which the distortion 
sets in and Ji + J2, one finds Tdist/{Ji + J2) = 0.5 ± 0.07 for Li2VOSi04 and a close value, 0.64± 0.07 for VOM0O4 
. This similarity can be understood by considering the expansion of the elastic and magnetic energies to lowest order 
in the displacements 

E {dJij/dXa){Si.Sj)Xa + ^ kafiXaXp/2 (9) 

with Xa.i3 the coordinates of the magnetic ions coupled by an elastic constant fc^jj. Since the reduction of magnetic 
energy is linear in Xa and the elastic one is quadratic, a minimum of magnetoelastic energy can be achieved for a 
small displacement Xeq of the coordinates. Now, if one considers just Ji and J2 couplings the order of magnitude 
of the energy gain induced by the displacement turns out to be E^^^ ~ —C[d{Ji + J2)/dXa]'^^^/kx^g, with C a 
constant which depends on the crystal structure. Then, if one considers that the similarities in Li2VOSi04 and 
VOM0O4 structure yield roughly similar elastic constants and power-law dependence of Ji on Xa, it is likely that 
ksTdist — Ex^g oc (Ji + J2), as experimentally found. 

It is remarkable to observe that while a clear signature of such a distortion is present in the NMR spectra, no 
modification in V'^'^ g-tensor is detected down to Tc- On the other hand, below 100 K a decrease in the intensity 
of the EPR signal, much faster than the decrease of the macroscopic magnetization, is observed (see Figs. || and ||). 
This effect is not associated with a saturation or a broadening of the EPR signal but rather indicates that there are 
some ions that are becoming EPR silent. These ions cannot correspond to V^^, which in a pyramidal coordination 
as the one in VOM0O4 , should always give an EPR signal. However, if valence fluctuations take place, they could 
correspond to Mo^"*" ions. In fact, due to selection rules, the signal of Mo^"*" in a regular tetrahedral configuration 
is cancelled out. The regular tetrahedral coordination is indeed supported by the small values of ^^Mo quadrupolar 
frequency (see Sect. IIB). As a whole, the comparison of NMR and EPR spectra leads to the following possible 
scenario. The distortion induced by the frustration causes, thanks to the hybridization of Mo d-orbitals in the band 
formation, a charge transfer from V^^ to Mo^+. As the distortion develops it induces a modification just in the 
NMR spectra of the adjacent nuclei yielding a broadening of the NMR line (see Sect. IIB) and the disappearence 
of the EPR signal of the adjacent ions. The V"^^ ions far from the distortion continue to give rise to an EPR 
signal with unchanged g-values, as experimentally observed. As the temperature is lowered the size of the distorted 
domains progressively grows and the EPR signal diminshes. The formation of distorted and non-distorted domains 
would support also the modifications in ^^Si NMR spectra in Li2VOSi04lj. In fact, in Li2VOSi04 as the temperature 
is lowered below Tdist one observes the progressive decrease of a low frequency peak (undistorted site) and the growth 
of a shifted high frequency peak (distorted site) (see Ref. ||). Hence, at Tdist a diffusive transition sets in yielding a 
progressive distortion of the whole lattice as the temperature decreases below 100 K. 

It should be noticed that it is somewhat unusual that Mo^"*" formation does not cause any Jahn- Teller effect. 
However, in this system the Jahn- Teller distortion could actually be hindered by the frustration driven distortion. 
This can occur if the distortion yields an energy gain of the frustrated magnetic lattice larger than the shift of the t2g 
ground-state. The charge transfer could also induce a progressive crossover of VOM0O4 from a half filled to a quarter 
filled band configuration, as the one of NaV2 05Ej, as the temperature decreases below 100 K and modifications in 
the transport properties should be observed. In fact, below 100 K a decrease in the energy barrier measured with 
resistivity is detectecCZl. 

IV. CONCLUSIONS 

In conlusion, it was shown that VOM0O4 is a prototype of a 2DFQHAF on a square lattice with Ji ~ J2, 
as Li2VOSi04 . The exchange couplings in VOM0O4 are much larger than the ones of Li2VOSi04 and a value 
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Ji + J2 — 155 K was derived, in good agreement with the one estimated from electronic structure calculations. In 
VOM0O4 a lattice distortion takes place at T^isf ~ 100 K. As the temperature is lowered below T^ist a progressive 
growth of the domains with lattice distortion occurs. From the comparison with Li2VOSi04 one finds that in these 
2DFQHAF T^ist roughly scales with Ji + J2, supporting the assumption that the distortion is driven by the magnetic 
frustration. Finally, in VOM0O4 novel phenomena, not observed in Li2VOSi04 , occur below Tdist and are tentatively 
associated with a charge transfer from V to Mo. 
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APPENDIX A: TIGHT BINDING MODEL 



The simplest model which can be used to describe the bands derived in the framework of the LDA includes 
twelve orbitals. The dxy of Vi and V2 (energy ed), where Vi is the atom at (1/4, 1/4, zy) and Vi is the atom at 
(-1/4, -1/4, -z^/); the 4 O p^ orbitals (energy Cp) centered in (±(1/4 - Ji), =F<52, T-^o), and (±(1/4 + =Fl/2 ± 
S2,Tzo)] the 4 O py orbitals (energy ep) centered in (=f5i, ±(1/4 + 52), T^^o) and ((=Fl/2 ± (5i), ±(1/4 — (52), T^o); 
the two Mo centered in (1/4, —1/4, 1/2) and (—1/4, 1/4, 1/2) (energy e,„). Within this model the following hopping 
integrals are considered: the hoppings between V d and O p states {tp^), the hopping between NN O p^ and O Py 
orbitals (too) and the hopping between Mo dxy and its NN O px and O py orbitals (tmo)- Starting from this bapd 
model the dispersion curve of the two conduction band can be obtained by downfolding all the O and V statesE3. 
Setting 5i = 62 = and neglecting t± one has 

e = £d + S-^ +4. ^^" ^_ (cos^(fcj2) + cos^(fc,/2)) 

^gjpd ^cos(fc,/2)cos(A:„/2), (Al) 

e — Cpl — 



with 



and 



b = 16 — "-^ \- 4 -1112. 1 + 4 — ::2£_ 

(e-ep)2 - ep){e - e,n) \ (e - Ep) 



'If / f 

^^no , mo ' - ■ 



1 + 4- 



(e-fp) (e - ep)(e - e„) V (e - Cp) 
Then one has that the in-plane NN and NNN hoppings are ti — 2^^^^ and ^2 — '^jzr^jzThi so that J2I J\ 
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FIG. 1: Structure of VOM0O4 projected along [001] . VO5 pyramids (black) run parallel to the line of sight and are connected 
by M0O4 tetrahedra (gray). For details see Ref. hd The dashed line shows the projection of the unit cell with a — 6.6078 A. 



FIG. 2: a) Derivative of the EPR powder spectra (solid lines) at T = 293 and 70 K. The circles show the best fit of the two 
spectra, used to determine g and AH . Notice that, in order to better illustrate the diferrent asymmetry, the intensity of two 
spectra was not normalized, b) Derivative of the EPR powder spectra for five selected temperatures between 69 K and Tc. The 
temperatures, moving from the most intense to the less intense signal, are T = 69, 59, 49, 45.5 and 42 K. A marked decrease in 
the intensity of the EPR signal on cooling is evident. 



FIG. 4: Temperature dependence of the area of the EPR signal (squares) and of the spin susceptibility (circles) measured 
experimentally with a SQUID magnetometer after subtraction of the Van-Vleck contribution (see text). The area of the EPR 
signal, which in principle is proportional to the spin susceptibility, was rescaled to match the value of the spin susceptibility at 
room temperature. 

FIG. 5: a) Temperature dependence of the susceptibility in VOM0O4 , for H — 1 kGauss. The solid line shows the high 
temperature Curie-Weiss behaviour for O = 155 K. b) Plot of the susceptibility measured with the SQUID magnetometer 
(after subtraction of the atomic diamagnetic contribution) versus the area of the EPR signal, with the temperature as an 
implicit parameter. The intercept of the solid line was used to derive Van-Vleck susceptibility. 



FIG. 6: a) ^^Mo NMR powder spectra in VOM0O4 oriented powders for = 9 Tesla along the c axis, b) ^^Mo NMR shift in 
VOM0O4 for H II c versus the spin susceptibility, measured with the SQUID magnetometer, after subtraction of the Van-Vleck 
term. The temperature, which is an implicit parameter, is shown for a few selected points. The solid lines evidence the change 
of slope, i.e. of hyperfine coupling, on cooling. 



FIG. 7: a) Temperature dependence of the resonance frequency for ®^Mo NMR central line, b) ®^Mo NMR shift of the central 
line in VOM0O4 unoriented powders plotted against the the spin susceptibility, measured with the SQUID magnetometer, after 
subtraction of the Van-Vleck term. The temperature, which is an implicit parameter, is shown for a few selected points. The 
solid lines evidence the change of slope, i.e. of hyperfine coupling, on cooling. 

FIG. 8: Temperature dependence of ^^Mo NMR l/Ti in VOM0O4 for the central line in a magnetic field of 9 Tesla, derived by 
fitting the recovery of nuclear magnetization with Eq. (^). 



FIG. 9: Temperature dependence of I/T2 in VOM0O4 powders for H = 9 Tesla, derived by fitting the echo decay of ®^Mo 
central transition with Eq. (m. 



FIG. 10: (Top) Band structure (left) and density of states (right) of VOM0O4 . The Fermi level is set at zero energy. The 
symmetry points are: r=(0,0,0), X=(7r/a,0,0), M=(7r/a, 7r/a,0), Z=(0,0,7r/c). (Bottom) Band structure of VOM0O4 close to 
the Fermi level. On the right side the total DOS (full line) is shown. The Vd projected DOS (dashed line), the Op and Mod 
projected DOS (dash-dotted and dotted lines) are also shown. 



FIG. 11: Spin susceptibility of Li2VOSi04 and VOM0O4 as a function of T/O, with 6 = 8.7 K and 155 K, respectively. The 
amplitude of the susceptibility of VOM0O4 has been rescaled by a factor slightly larger than the ratio between the Curie- Weiss 
temperatures indicating a slightly lower purity of VOM0O4 sample with respect to Li2VOSi04 . 

FIG. 12: Temperature dependence of ^^V NMR l/Ti as estimated from the temperature dependence of ^®Mo I/T2 shown in 
Fig. ^ (see Sect. IIIB). The line is a guide to the eye. 



FIG. 3: Temperature dependence of the peak to peak width of the EPR powder spectrum in VOM0O4 . In the inset the ratio 
between the EPR linewidth for H \\ and _L c, estimated from the analysis of the powder spectra, is reported. 
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